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Abstract. 

Relatively uniform light curves and spectral evolution of Type la supernovae (SNe 
la) have led to the use of SNe la as a "standard candle" to determine cosmological 
parameters. Whether a statistically significant value of the cosmological constant can 
be obtained depends on whether the peak luminosities of SNe la are sufficiently free 
from the effects of cosmic and galactic evolutions. 

Here we first review the single degenerate scenario for the Chandrasekhar mass white 
dwarf (WD) models of SNe la. We identify the progenitor's evolution and population 
with two channels: (1) the WD+RG (red-giant) and (2) the WD+MS (near main- 
sequence He-rich star) channels. In these channels, the strong wind from accreting WDs 
plays a key role, which yields important age and metallicity effects on the evolution. 

We then address the questions whether the nature of SNe la depends systematically 
on environmental properties such as metallicity and age of the progenitor system and 
whether significant evolutionary effects exist. We suggest that the variation of the 
carbon mass fraction X(C) in the C+O WD (or the variation of the initial WD mass) 
causes the diversity of the brightness of SNe la. This model can explain the observed 
dependences of SNe la brightness on the galaxy types and the distance from the galactic 
center. 

Finally, applying the metallicity effect on the evolution of SN la progenitors, we make 
a prediction of the cosmic supernova rate history as a composite of the supernova rates 
in different types of galaxies. 



I INTRODUCTION 

Type la supernovae (SNe la) are good distance indicators, and provide a promis- 
ing tool for determining cosmological parameters (e.g., [1]). SNe la have been 
discovered up to z ~ 1.32 [2]. Both the Supernova Cosmology Project [3,4] and the 
High-z Supernova Search Team [5,6] have suggested a statistically significant value 
for the cosmological constant. 



However, SNe la are not perfect standard candles, but show some intrinsic varia- 
tions in brightness. When determining the absolute peak luminosity of high-redshift 
SNe la, therefore, these analyses have taken advantage of the empirical relation ex- 
isting between the peak brightness and the light curve shape (LCS). Since this re- 
lation has been obtained from nearby SNe la only [7-9], it is important to examine 
whether it depends systematically on environmental properties such as metallicity 
and age of the progenitor system. 

High-redshift supernovae present us very useful information, not only to deter- 
mine cosmological parameters but also to put constraints on the star formation 
history in the universe. They have given the SN la rate at z ~ 0.5 [10] but will 
provide the SN la rate history over < z < 1. With the Next Generation Space 
Telescope, both SNe la and SNe II will be observed through z ~ 4. It is useful to 
provide a prediction of cosmic supernova rates to constrain the age and metallicity 
effects of the SN la progenitors. 

SNe la have been widely believed to be a thermonuclear explosion of a mass- 
accreting white dwarf (WD) (e.g., [11] for a review). However, the immediate 
progenitor binary systems have not been clearly identified yet [12]. In order to 
address the above questions regarding the nature of high-redshift SNe la, we need 
to identify the progenitors systems and examine the "evolutionary" effects (or en- 
vironmental effects) on those systems. 

In §2, we summarize the progenitors' evolution where the strong wind from ac- 
creting WDs plays a key role [13-15]. In §3, we address the issue of whether a 
difference in the environmental properties is at the basis of the observed range of 
peak brightness [17]. In §4, we make a prediction of the cosmic supernova rate 
history as a composite of the different types of galaxies [18]. 

II EVOLUTION OF PROGENITOR SYSTEMS 

There exist two models proposed as progenitors of SNe la: 1) the Chandrasekhar 
mass model, in which a mass-accreting carbon-oxygen (C+O) WD grows in mass 
up to the critical mass Mj a ~ 1.37 — 1.38M Q near the Chandrasekhar mass and 
explodes as an SN la (e.g., [19,20]), and 2) the sub- Chandrasekhar mass model, 
in which an accreted layer of helium atop a C+O WD ignites off-center for a WD 
mass well below the Chandrasekhar mass (e.g., [21]). The early time spectra of 
the majority of SNe la are in excellent agreement with the synthetic spectra of 
the Chandrasekhar mass models, while the spectra of the sub- Chandrasekhar mass 
models are too blue to be comparable with the observations [22,23]. 

For the evolution of accreting WDs toward the Chandrasekhar mass, two sce- 
narios have been proposed: 1) a double degenerate (DD) scenario, i.e., merging of 
double C+O WDs with a combined mass surpassing the Chandrasekhar mass limit 
[24,25], and 2) a single degenerate (SD) scenario, i.e., accretion of hydrogen-rich 
matter via mass transfer from a binary companion (e.g., [26,20]). The issue of DD 
vs. SD is still debated (e.g., [12]), although theoretical modeling has indicated 



that the merging of WDs leads to the accretion-induced collapse rather than SN la 
explosion [27-29]. 

In the SD Chandrasekhar mass model for SNe la, a WD explodes as a SN la only 
when its rate of the mass accretion (M) is in a certain narrow range (e.g., [26,30]). 
In particular, if M exceeds the critical rate M b , the accreted matter extends to 
form a common envelope [31]. This difficulty has been overcome by the WD wind 
model (see below). For the actual binary systems which grow the WD mass (Mwd) 
to Mj a , the following two systems are appropriate. One is a system consisting of a 
mass-accreting WD and a lobe-filling, more massive, slightly evolved main-sequence 
or sub-giant star (hereafter "WD+MS system"). The other system consists of a 
WD and a lobe-filling, less massive, red-giant (hereafter "WD+RG system"). 

A White dwarf winds 

Optically thick WD winds are driven when the accretion rate M exceeds the 
critical rate M b . Here M b is the rate at which steady burning can process the 
accreted hydrogen into helium as M b « 0.75 x 10~ 6 (jjgg- - 0.40) M yr" 1 . 

With such a rapid accretion, the WD envelope expands to i? p h ~ 0.1 R@ and the 
photospheric temperature decreases below logT p h ~ 5.5. Around this temperature, 
the shoulder of the strong peak of OPAL Fe opacity [32] drives the radiation- 
driven wind [13,15]. The ratio of f p h/f C sc between the photospheric velocity and 
the escape velocity at the photosphere depends on the mass transfer rate and Mwd- 
(see Fig. 6 in [15]). We call the wind strong when v p h > v esc . When the wind is 
strong, f ph ~ 1000 km s _1 being much faster than the orbital velocity. 

If the wind is sufficiently strong, the WD can avoid the formation of a common 
envelope and steady hydrogen burning increases its mass continuously at a rate M b 
by blowing the extra mass away in a wind. When the mass transfer rate decreases 
below this critical value, optically thick winds stop. If the mass transfer rate further 
decreases below ~ 0.5 M b , hydrogen shell burning becomes unstable to trigger very 
weak shell flashes but still burns a large fraction of accreted hydrogen. 

The steady hydrogen shell burning converts hydrogen into helium atop the C+O 
core and increases the mass of the helium layer gradually. When its mass reaches a 
certain value, weak helium shell flashes occur. Then a part of the envelope mass is 
blown off but a large fraction of He can be burned to C+O [33] to increase the WD 
mass. In this way, strong winds from the accreting WD play a key role to increase 
the WD mass to M Ia . 

B WD+RG system 

This is a symbiotic binary system consisting of a WD and a low mass red-giant 
(RG). A full evolutionary path of the WD+RG system from the zero age main- 
sequence stage to the SN la explosion is described in [15,16]. The occurrence 
frequency of SNe la through this channel is much larger than the earlier scenario, 
because of the following two evolutionary processes, which have not considered 
before. 



(1) Because of the AGB wind, the WD + RG close binary can form from a wide 
binary even with such a large initial separation as ai < 40, OOOi? . Our earlier 
estimate [13] is constrained by a, < l,5OOi? . 

(2) When the RG fills its inner critical Roche lobe, the WD undergoes rapid 
mass accretion and blows a strong optically thick wind. Our earlier analysis has 
shown that the mass transfer is stabilized by this wind only when the mass ratio 
of RG/WD is smaller than 1.15. Our new finding is that the WD wind can strip 
mass from the RG envelope, which could be efficient enough to stabilize the mass 
transfer even if the RG/WD mass ratio exceeds 1.15. If this mass-stripping effect 
is strong enough, though its efficiency i] cS is subject to uncertainties, the symbiotic 
channel can produce SNe la for a much (ten times or more) wider range of the 
binary parameters than our earlier estimation. 

With the above two new effects (1) and (2), the WD+RG (symbiotic) channel 
can account for the inferred rate of SNe la in our Galaxy. The immediate pro- 
genitor binaries in this symbiotic channel to SNe la may be observed as symbiotic 
stars, luminous supersoft X-ray sources, or recurrent novae like T CrB or RS Oph, 
depending on the wind status. 

C WD+MS system 

In this scenario, a C+O WD is originated, not from an AGB star with a C+O 
core, but from a red-giant star with a helium core of ~ 0.8 — 2.OM . The helium 
star, which is formed after the first common envelope evolution, evolves to form a 
C+O WD of ~ 0.8 — I.IMq with transferring a part of the helium envelope onto 
the secondary main-sequence star. A full evolutionary path of the WD+MS system 
from the zero age main-sequence stage to the SN la explosion is described in [14]. 

This evolutionary path provides a much wider channel to SNe la than previous 
scenarios. A part of the progenitor systems are identified as the luminous supersoft 
X-ray sources [34] during steady H-burning (but without wind to avoid extinction), 
or the recurrent novae like U Sco if H-burning is weakly unstable. Actually these 
objects are characterized by the accretion of helium-rich matter. 

D Realization frequency 

For an immediate progenitor system WD+RG of SNe la, we consider a close 
binary initially consisting of a C+O WD with Mwd.o = 0.6 — 1.2M and a low- 
mass red-giant star with M RG ,o = 0.7 — 3.OM having a helium core of M Hc ,o = 
0.2 — O.46M . The initial state of these immediate progenitors is specified by three 
parameters, i.e., Mwd.o, ^rg,o = M^o, and the initial orbital period Po (MH e ,o is 
determined if P is given). 

We follow binary evolutions of these systems and obtain the parameter range(s) 
which can produce an SN la. In Figure 1, the region enclosed by the thin solid line 
produces SNe la for several cases of the initial WD mass, Mwd,o — 0.75 - 1.1 M & . 
For smaller M W d,o> the wind is weaker, so that the SN la region is smaller. The 
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FIGURE 1. The region to produce SNe la in the logPo — M^o plane for five initial WD masses 
of O.75M , O.8M , O.9M , l.OM (heavy solid line), and 1.1M Q . The region of M W d,o = O.7M 
almost vanishes for both the WD+MS and WD+RG systems, and the region of Mwd,o = O.75M 
vanishes for the WD+RG system. Here, we assume the stripping efficiency of ?7 c ff = 1. For 
comparison, we show only the region of -Mwd.o = 1-OM for a much lower efficiency of ?7 e ff = 0.3 
by a dash-dotted line. 

regions of M WD ,o = O.6M and O.7M vanish for both the WD+MS and WD+RG 
systems. 

In the outside of this region, the outcome of the evolution at the end of the 
calculations is not an SN la but one of the followings: (i) Formation of a common 
envelope for too large M d or P ~ day, where the mass transfer is unstable at the 
beginning of mass transfer, (ii) Novae or strong hydrogen shell flash for too small 
Md,o, where the mass transfer rate becomes below KT 7 M Q yr l . (iii) Helium core 
flash of the red giant component for too long P , where a central helium core flash 
ignites, i.e., the helium core mass of the red-giant reaches O.46M . (iv) Accretion- 
induced collapse for Mwd,o > 1-2M , where the central density of the WD reaches 
~ 10 10 g cm -3 before heating wave from the hydrogen burning layer reaches the 
center. As a result, the WD undergoes collapse due to electron capture without 
exploding as an SN la [30]. 

It is clear that the new region of the WD+RG system is not limited by the 
condition of q < 1.15, thus being ten times or more wider than the region of [13] 's 
model (depending on the the stripping efficiency of r] c g). 

The WD+MS progenitor system can also be specified by three initial parameters: 
the initial C+O WD mass M WDj0 , the mass donor's initial mass M dj0 , and the 
orbital period P . For M W d,o — 1-OM , the region producing an SN la is bounded 
by M&fl = 1.8 — 3.2M and Pq = 0.5 — 5 d as shown by the solid line in Figure 1. 
The upper and lower bounds are respectively determined by the common envelope 
formation (i) and nova-like explosions (ii) as above. The left and right bounds are 
determined by the minimum and maximum radii during the main sequence of the 
donor star [14]. 

We estimate the rate of SNe la originating from these channels in our Galaxy 
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FIGURE 2. The regions of SNc la is plot- 
ted in the initial orbital period vs. the initial 
companion mass diagram for the initial WD 
mass of M\vd,o = 1.0M Q . The dashed and 
solid lines represent the cases of solar abun- 
dance (Z = 0.02) and much lower metallicity 
of Z = 0.004, respectively. The left and the 
right regions correspond to the WD+MS and 
the WD+RG systems, respectively. 
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FIGURE 3. WD mass vs. metallicity dia- 
gram showing the metallicity dependence of 
optically thick winds. We regard the wind as 
"strong" if v w > v csc but "weak" if v w < v csc . 
The term of "He" or "Fe" wind denotes that 
the wind is accelerated by the peak of iron 
lines near logT(K) ~ 5.2 or of helium lines 
near logT(K) ~ 4.6. The dashed line indi- 
cates the demarcation between the "strong" 
wind and the "weak" wind. 



by using equation (1) of [24]. The realization frequencies of SNe la through the 
WD+RG and WD+MS channels are estimated as ~ 0.0017 yr" 1 (WD+RG) and ~ 
0.001 yr" 1 (WD+MS), respectively. The total SN la rate of the WD+MS/WD+RG 
systems becomes ~ 0.003 yr~ x , which is close enough to the inferred rate of our 
Galaxy. 

E Low metallicity inhibition of type la supernovae 

The optically thick winds are driven by a strong peak of OPAL opacity at 
logT(K) ~ 5.2 (e.g., [32]). Since the opacity peak is due to iron lines, the wind 
velocity v w depends on the iron abundance [Fe/H] ( [35,36]), i.e., v w is higher for 
larger [Fe/H]. The metallicity effect on SNe la is clearly demonstrated by the size 
of the regions to produce SNe la in the diagram of the initial orbital period versus 
initial mass of the companion star (see Fig. 2). The SN la regions are much smaller 
for lower metallicity because the wind becomes weaker. 

The wind velocity depends also on the luminosity L of the WD. The more massive 
WD has a higher L, thus blowing higher velocity winds ( [15]). In order for the 
wind velocity to exceed the escape velocity of the WD near the photosphere, the 
WD mass should be larger than a certain critical mass for a given [Fe/H]. This 
implies that the initial mass of the WD M WDj0 should already exceed that critical 
mass in order for the WD mass to grow to the Ch mass. This critical mass is larger 
for smaller [Fe/H], reaching 1.1M for [Fe/H] = —1.1 (Fig. 3). Here we should 
note that the relative number of WDs with M WD > 1.1 M© is quite small in close 



binary systems ( [47]). And for M WDi0 > 1.2M , the accretion leads to collapse 
rather than SNe la ( [30]). Therefore, no SN la occurs at [Fe/H] < —1.1 in our 
model. 

It is possible to test the metallicity effects on SNe la with the chemical evolution 
of galaxies. 

In the one-zone uniform model for the chemical evolution of the solar neighbor- 
hood, the heavy elements in the metal-poor stars originate from the mixture of 
the SN II ejecta of various progenitor masses. The abundances averaged over the 
progenitor masses of SNe II predicts [O/Fe] ~ 0.45 (e.g., [37,38]). Later SNe la 
start ejecting mostly Fe, so that [O/Fe] decreases to ~ around [Fe/H] ~ 0. The 
low-metallicity inhibition of SNe la predicts that the decrease in [O/Fe] starts at 
[Fe/H] ~ —1. Such an evolution of [O/Fe] well explains the observations ( [35]). 

However, we should note that some anomalous stars have [O/Fe] ~ at [Fe/H] 
< — 1. The presence of such stars, however, is not in conflict with our SNe la 
models, but can be understood as follows: The formation of such anomalous stars 
(and the diversity of [O/Fe] in general) indicates that the interstellar materials were 
not uniformly mixed but contaminated by only a few SNe II (or even single SN II) 
ejecta. This is because the timescale of mixing was longer than the time difference 
between the supernova event and the next generation star formation. The iron 
and oxygen abundances produced by a single SN II vary depending on the mass, 
energy, mass cut, and metallicity of the progenitor. Relatively smaller mass SNe II 
(13 — 15M Q ) and higher explosion energies tend to produce [O/Fe] ~ ( [38,39]). 
Those metal poor stars with [O /Fe] ~ may be born from the interstellar medium 
polluted by such SNe II. 

The metallicity effect on SNe la can also be checked with the metallicity of the 
host galaxies of nearby SNe la. There has been no evidence that SNe la have 
occurred in galaxies with a metallicity of [Fe/H] < —1, although host galaxies 
are detected only for one third of SNe la and the estimated metallicities of host 
galaxies are uncertain. Three SNe la are observed in low-metallicity dwarf galaxies; 
SN1895B and SN1972E in NGC 5253, and SN1937C in IC 4182. Metallicities of 
these galaxies are estimated to be [O/H] = —0.25 and —0.35, respectively [40]. If 
[O/Fe] ~ as in the Magellanic Clouds, [Fe/H]~ —0.25 and —0.35 which are not 
so small. Even if these galaxies have extremely SN II like abundance as [O/Fe] 
~ 0.45, [Fe/H] ~ —0.7 and —0.8 (being higher than —1), respectively. Since these 
host galaxies are blue (B — V — 0.44 for NGC 5253 and B — V = 0.37 for IC 4182 
according to RC3 catalog), the MS+WD systems are dominant progenitors for the 
present SNe la. The rate of SNe la originated from the MS+WD systems is not so 
sensitive to the metallicity as far as [Fe/H] > —1 ( [36]). Even if [Fe/H] ~ —0.7 in 
such blue galaxies, therefore, the SN la rate is predicted to be similar to those in 
more metal-rich galaxies. 
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FIGURE 4. Abundances in mass fraction in the inner core of the 6 Mq star for Y = 0.2775 
and Z = 0.02 at the end of the second dredge-up. 

Ill THE ORIGIN OF DIVERSITY OF SNE IA AND 
ENVIRONMENTAL EFFECTS 

There are some observational indications that SNe la are affected by their envi- 
ronment. The most luminous SNe la seem to occur only in spiral galaxies, while 
both spiral and elliptical galaxies are hosts for dimmer SNe la. Thus the mean 
peak brightness is dimmer in ellipticals than in spiral galaxies [41]. The SNe la 
rate per unit luminosity at the present epoch is almost twice as high in spirals as 
in ellipticals [42]. Moreover, [43,44] found that the variation of the peak brightness 
for SNe located in the outer regions in galaxies is smaller. 

[45,46] examined how the initial composition of the WD (metallicity and the 
C/O ratio) affects the observed properties of SNe la. [47] obtained the C/O ratio 
as a function of the main-sequence mass and metallicity of the WD progenitors. 
[17] suggested that the variation of the C/O ratio is the main cause of the variation 
of SNe la brightness, with larger C/O ratio yielding brighter SNe la. We will 
show that the C/O ratio depends indeed on environmental properties, such as the 
metallicity and age of the companion of the WD, and that our model can explain 
most of the observational trends discussed above. We then make some predictions 
about the brightness of SN la at higher redshift. 

A C/O ratio in WD progenitors 

In this section we discuss how the C/O ratio in the WD depends on the metallicity 
and age of the binary system. The C/O ratio in C+O WDs depends primarily on 
the main-sequence mass of the WD progenitor and on metallicity. 



We calculated the evolution of intermediate- mass (3 — 9M ) stars for metallicity 
Z=0.001 - 0.03. In the ranges of stellar masses and Z considered in this paper, 
the most important metallicity effect is that the radiative opacity is smaller for 
lower Z . Therefore, a star with lower Z is brighter, thus having a shorter lifetime 
than a star with the same mass but higher Z . In this sense, the effect of reducing 
metallicity for these stars is almost equivalent to increasing a stellar mass. 

For stars with larger masses and/or smaller Z, the luminosity is higher at the 
same evolutionary phase. With a higher nuclear energy generation rate, these stars 
have larger convective cores during H and He burning, thus forming larger He and 
C-0 cores. 

As seen in Figure 4, the central part of these stars is oxygen-rich. The C/O ratio 
is nearly constant in the innermost region, which was a convective core during He 
burning. Outside this homogeneous region, where the C-0 layer grows due to He 
shell burning, the C/O ratio increases up to C/O > 1; thus the oxygen-rich core is 
surrounded by a shell with C/O > 1. In fact this is a generic feature in all models 
we calculated. The C/O ratio in the shell is C/O ~ 1 for the star as massive as 
~ 7M Q , and C/O > 1 for less massive stars. 

When a progenitor reaches the critical mass for the SNe la explosion, the central 
core is convective up to around 1.1 M Q . Hence the relevant C/O ratio is between 
the central value before convective mixing and the total C/O of the whole WD. 
Using the results from the C6 model [19], we assume that the convective region is 
1.14 M and for simplicity, C/O = 1 outside the C-0 core at the end of second 
dredge-up. Then we obtain the C/O ratio of the inner part of the SNe la progenitors 
(Fig- 5). 

From this figure we find three interesting trends. First, while the central C/O is 
a complicated function of stellar mass [47], as shown here the C/O ratio in the core 
before SNe la explosion is a decreasing monotonic function of mass. The central 
C/O ratio at the end of second dredge-up decreases with mass for M ms > 5M , 
while the ratio increases with mass for M ms > 4M Q ; however, the convective core 
mass during He burning is smaller for a less massive star, and the C/O ratio during 
shell He burning is larger for smaller C+O core. Hence, when the C/O ratio is 
averaged over 1.1 M Q the C/O ratio decreases with mass. Second, as shown in 
[47], although the C/O ratio is a complicated function of metallicity and mass, the 
metallicity dependence is remarkably converged when the ratio is seen as a function 
of the C+O core mass (Mqo) instead of the initial main sequence mass. 

According to the evolutionary calculations for 3—9 M stars by [47], the C/O 
ratio and its distribution are determined in the following evolutionary stages of the 
close binary. 

(1) At the end of central He burning in the 3—9 M primary star, C/0< 1 in 
the convective core. The mass of the core is larger for more massive stars. 

(2) After central He exhaustion, the outer C+O layer grows via He shell burning, 
where C/0> 1 [47]. 

(3a) If the primary star becomes a red giant (case C evolution; e.g., [48]), it 
then undergoes the second dredge-up, forming a thin He layer, and enters the AGB 
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FIGURE 5. The total 12 C mass fraction included in the convective core of mass, M = 1.14M , 
just before the SN la explosion as a function of the C+O core mass before the onset of mass 
accretion, Mco- The lower bounds of Mqo obtained from the age effects and the conditions for 
strong wind to blow are also shown by arrows. 

phase. The C+O core mass, M co , at this phase is larger for more massive stars. 
For a larger M co the total carbon mass fraction is smaller. 

(3b) When it enters the AGB phase, the star greatly expands and is assumed 
here to undergo Roche lobe overflow (or a super-wind phase) and to form a C+O 
WD. Thus the initial mass of the WD, M W d,o> in the close binary at the beginning 
of mass accretion is approximately equal to M co . 

(4a) If the primary star becomes a He star (case BB evolution), the second 
dredge-up in (3a) corresponds to the expansion of the He envelope. 

(4b) The ensuing Roche lobe overflow again leads to a WD of mass Mwd,o = 
M co . 

(5) After the onset of mass accretion, the WD mass grows through steady H 
burning and weak He shell flashes, as described in the WD wind model. The 
composition of the growing C+O layer is assumed to be C/0=1. 

(6) The WD grows in mass and ignites carbon when its mass reaches Mi a = 
1.367M Q , as in the model C6 of [19]. Because of strong electron-degeneracy carbon 
burning is unstable and grows into a deflagration for a central temperature of 8 x 10 8 
K and a central density of 1.47 x 10 9 g cm -3 . At this stage, the convective core 
extends to M r = 1.14M and the material is mixed almost uniformly as in the C6 
model. 

In Figure 5, we show the carbon mass fraction X(C) in the convective core of 
this pre-explosive WD, as a function of metallicity (Z) and initial mass of the WD 
before the onset of mass accretion, Mqo- Figure 5 reveals that: 1) X(C) is smaller 
for larger Mqo — ^wd,o- 2) The dependence of X(C) on metallicity is small when 
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plotted against Mco, even though the relation between Mqo and the initial stellar 
mass depends sensitively on Z [47]. 



B Brightness of SNe la and the C/O ratio 

In the Chandrasekhar mass models for SNe la, the brightness of SNe la is deter- 
mined mainly by the mass of 56 Ni synthesized (M^q). Observational data suggest 
that M Ni56 for most SNe la lies in the range M Nm ~ 0.4 - O.8M (e.g., [49]). This 
range of Mjsfi56 can result from differences in the C/O ratio in the progenitor WD 
as follows. 

In the deflagration model, a larger C/O ratio leads to the production of more 
nuclear energy and buoyancy force, thus leading to a faster propagation. The 
faster propagation of the convective deflagration wave results in a larger M^ 5e . 
For example, a variation of the propagation speed by 15% in the W6 - W8 models 
results in M Ni56 values ranging between 0.5 and O.7M [19], which could explain 
the observations. 

In the delayed detonation model, is predominantly determined by the 

deflagration-to-detonation-transition (DDT) density pddt, at which the initially 
subsonic deflagration turns into a supersonic detonation [50]. As discussed in [17], 
Pddt could be very sensitive to X(C), and a larger X(C) is likely to result in a 
larger p DDT and M Ni56 . 

Here we postulate that Mn;56 and consequently brightness of a SN la increase as 
the progenitors' C/O ratio increases (and thus Mwd,o decreases). As illustrated in 
Figure 5, the range of M Ni56 ~ 0.5 — O.8M is the result of an X(C) range 0.35 — 0.5, 
which is the range of X(C) values of our progenitor models. The X(C) - M Ni56 
- M WD relation we adopt is still only a working hypothesis, which needs to be 
proved from studies of the turbulent flame during explosion (e.g., [51]). 



C Metallicity and age effects 

1 Metallicity effects on the minimum Mwd,o 

As mentioned in §2.5, M w is the metallicity-dependent minimum M W d,o f° r a 
WD to become an SN la (strong wind condition in Fig. 5). The upper bound 
M WDj0 ~ l.O7M is imposed by the condition that carbon should not ignite and is 
almost independent of metallicity. As shown in Figure 5, the range of Mqq ~ Mwd,o 
can be converted into a range of X(C). From this we find the following metallicity 
dependence for X(C): 

(1) The upper bound of X(C), which is determined by the lower limit on Mqo 
imposed by the metallicity-dependent conditions for a strong wind, e.g., X(C) 
< 0.51, 0.46 and 0.41, for Z=0.02, 0.01, and 0.004, respectively. 

(2) On the other hand, the lower bound, X(C) ~ 0.35 — 0.33, does not depend 
much on Z, since it is imposed by the maximum M co . 

(3) Assuming the relation between M Ni56 and X(C) given in Figure 5, our model 
predicts the absence of brighter SNe la in lower metallicity environment. 

2 Age effects on the minimum Mwd,o 

In our model, the age of the progenitor system also constrains the range of X(C) 
in SNe la. In the SD scenario, the lifetime of the binary system is essentially the 
main-sequence lifetime of the companion star, which depends on its initial mass M 2 . 
[14,15] have obtained a constraint on M 2 by calculating the evolution of accreting 
WDs for a set of initial masses of the WD (Mwd,o — Mco) and of the companion 
(M2), and the initial binary period (Pq). In order for the WD mass to reach Mi a , 
the donor star should transfer enough material at the appropriate accretion rates. 
The donors of successful cases are divided into two categories: one is composed of 
slightly evolved main-sequence stars with M 2 ~ 1.7 — 3.6M Q (for Z=0.02), and the 
other of red-giant stars with M 2 ~ 0.8 - 3.1M© (for Z=0.02) (Fig. 1). 

If the progenitor system is older than 2 Gyr, it should be a system with a donor 
star of M 2 < 1.7 M® in the red-giant branch. Systems with M 2 > 1.7M Q become 
SNe la in a time shorter than 2 Gyr. Likewise, for a given age of the progenitor 
system, M 2 must be smaller than a limiting mass. This constraint on M 2 can be 
translated into the presence of a minimum Mco f° r a given age, as follows: For 
a smaller M 2 , i.e. for the older system, the total mass which can be transferred 
from the donor to the WD is smaller. In order for M WD to reach M Ia , therefore, 
the initial mass of the WD, M WD ~ M co , should be larger. This implies that 
the older system should have larger minimum Mqo as indicated in Figure 5. Using 
the X(C)-Mco and M-^ i5 Q-X(C) relations (Fig. 5), we conclude that WDs in older 
progenitor systems have a smaller X(C), and thus produce dimmer SNe la. 

D Comparison with observations 

The first observational indication which can be compared with our model is the 
possible dependence of the SN brightness on the morphology of the host galaxies. 



[41] found that the most luminous SNe la occur in spiral galaxies, while both spiral 
and elliptical galaxies are hosts to dimmer SNe la. Hence, the mean peak brightness 
is lower in elliptical than in spiral galaxies. 

In our model, this property is simply understood as the effect of the different age 
of the companion. In spiral galaxies, star formation occurs continuously up to the 
present time. Hence, both WD+MS and WD+RG systems can produce SNe la. 
In elliptical galaxies, on the other hand, star formation has long ended, typically 
more than 10 Gyr ago. Hence, WD+MS systems can no longer produce SNe la. 
In Figure 6, we show the frequency of the expected SN I for a galaxy of mass 
2 x 1O U M for WD+MS and WD+RG systems separately as a function of M co . 
Here we use the results of [15,14], and the M co —X(C) and M Ni56 — X(C) relations 
given in Figure 5. Since a WD with smaller Mqo is assumed to produce a brighter 
SN la (larger M^q), our model predicts that dimmer SNe la occur both in spirals 
and in ellipticals, while brighter ones occur only in spirals. The mean brightness 
is smaller for ellipticals and the total SN la rate per unit luminosity is larger in 
spirals than in ellipticals. These properties are consistent with observations. 

The second observational suggestion is the radial distribution of SNe la in galax- 
ies. [43,6] found that the variation of the peak brightness for SNe la located in the 
outer regions in galaxies is smaller. This behavior can be understood as the effect 
of metallicity. As shown in Figure 5, even when the progenitor age is the same, 
the minimum Mco is larger for a smaller metallicity because of the metallicity 
dependence of the WD winds. Therefore, our model predicts that the maximum 
brightness of SNe la decreases as metallicity decreases. Since the outer regions 
of galaxies are thought to have lower metallicities than the inner regions [52,53], 
our model is consistent with observations. [43] also claimed that SNe la may be 
deficient in the bulges of spiral galaxies. This can be explained by the age effect, 
because the bulge consists of old population stars. 

E Evolution of SNe la at high redshift 

We have suggested that X(C) is the quantity very likely to cause the diversity 
in M Ni56 and thus in the brightness of SNe la. We have then shown that our model 
predicts that the brightness of SNe la depends on the environment, in a way which 
is qualitatively consistent with the observations. Further studies of the propagation 
of the turbulent flame and the DDT are necessary in order to actually prove that 
X(C) is the key parameter. 

Our model predicts that when the progenitors belong to an old population, or 
to a low metal environment, the number of very bright SNe la is small, so that 
the variation in brightness is also smaller, which is shown in Figure 7. In spiral 
galaxies, the metallicity is significantly smaller at redshifts z J> 1, and thus both 
the mean brightness of SNe la and its range tend to be smaller (Fig. 7). At z > 2 
SNe la would not occur in spirals at all because the metallicity is too low. In 
elliptical galaxies, on the other hand, the metallicity at redshifts z ~ 1 — 3 is not 
very different from the present value. However, the age of the galaxies at z ~ 1 is 
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FIGURE 7. Illustration of the predicted variation in SN la brightness with redshift. 



only about 5 Gyr, so that the mean brightness of SNe la and its range tend to be 
larger at z > 1 than in the present ellipticals because of the age effect. 

We note that the variation of X(C) is larger in metal-rich nearby spirals than 
in high redshift galaxies. Therefore, if X(C) is the main parameter responsible for 
the diversity of SNe la, and if the LCS method is confirmed by the nearby SNe la 
data, the LCS method can also be used to determine the absolute magnitude of 
high redshift SNe la. 

F Possible evolutionary effects 

In the above subsections, we consider the metallicity effects only on the C/O 
ratio; this is just to shift the main-sequence mass - Mwd.o relation, thus resulting 
in no important evolutionary effect. However, some other metallicity effects could 
give rise to evolution of SNe la between high and low redshifts (i.e., between low 
and high metallicities). 

Here we point out just one possible metallicity effect on the carbon ignition 
density in the accreting WD. The ignition density is determined by the competition 
between the compressional heating due to accretion and the neutrino cooling. The 
neutrino emission is enhanced by the local Urea shell process of, e.g., 21 Ne- 21 F pair 
[54]. (Note that this is different from the convective Urea neutrino process). For 
higher metallicity, the abundance of 21 Ne is larger so that the cooling is larger. 
This could delay the carbon ignition until a higher central density is reached [55]. 

Since the WD with a higher central density has a larger binding energy, the 
kinetic energy of SNe la tends to be smaller if the same amount of 56 Ni is produced. 
This might cause a systematically slower light curve evolution at higher metallicity 
environment. The carbon ignition process including these metallicity effects as 
well as the convective Urea neutrino process need to be studied (see also [56] for 
nucleosynthesis constraints on the ignition density). 

IV COSMIC SUPERNOVA RATES 



Attempts have been made to predict the cosmic supernova rates as a function 
of redshift by using the observed cosmic star formation rate (SFR) [57-59]. The 
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observed cosmic SFR shows a peak at z ~ 1.4 and a sharp decrease to the present 
[60] . However, UV luminosities which is converted to the SFRs may be affected by 
the dust extinction [61]. Recent updates of the cosmic SFR suggest that a peak 
lies around z ~ 3. 

[35] predicts that the cosmic SN la rate drops at z ~ 1 — 2, due to the metallicity- 
dependence of the SN la rate. Their finding that the occurrence of SNe la depends 
on the metallicity of the progenitor systems implies that the SN la rate strongly 
depends on the history of the star formation and metal-enrichment. The universe 
is composed of different morphological types of galaxies and therefore the cosmic 
SFR is a sum of the SFRs for different types of galaxies. As each morphological 
type has a unique star formation history, we should decompose the cosmic SFR 
into the SFR belonging to each type of galaxy and calculate the SN la rate for each 
type of galaxy. 

Here we first construct the detailed evolution models for different type of galaxies 
which are compatible with the stringent observational constraints, and apply them 
to reproduce the cosmic SFR for two different environments, e.g., the cluster and 
the field. Secondly with the self-consistent galaxy models, we calculate the SN rate 
history for each type of galaxy and predict the cosmic supernova rates as a function 
of redshift. 



Galaxies that are responsible for the cosmic SFR have different timescales for the 
heavy-element enrichment, and the occurrence of supernovae depends on the metal- 
licity therein. Therefore we calculate the cosmic supernova rate by summing up 
the supernova rates in spirals (SOa-Sa, Sab-Sb, Sbc-Sc, and Scd-Sd) and ellipticals 
with the ratio of the relative mass contribution. The relative mass contribution is 



A In Clusters 



obtained from the observed relative luminosity proportion and the calculated mass 
to light ratio in B-band [18]. The photometric evolution is calculated with the 
spectral synthesis population database taken from [62]. We adopt H = 50 km s _1 
Mpc -1 , Qq = 0.2, Ao = 0, and the galactic age of t age = 15 Gyr. 

First, we make a prediction of the cosmic supernova rates by using the galaxy 
models which are constructed to meet the observational constraints of cluster galax- 
ies. We assume that elliptical galaxies are formed by a single star burst and stop 
the star formation at t ~ 1 Gyr due to a supernova-driven galactic wind, while spi- 
ral galaxies are formed by a relatively continuous star formation. The infall rates 
and the SFRs are given by [18]. These models are constructed to meet the latest 
observational constraints such as the present gas fractions and colors for spirals, 
and the mean stellar metallicity and the color evolution from the present to z ~ 1 
for ellipticals (see [18] for the figures). 

The synthesized cosmic SFR has an excess at z > 3 due to the early star burst in 
ellipticals and a shallower slope from the present to the peak at z ~ 1.4, compared 
with Madau's plot [60]. Figure 8 shows the cosmic supernova rates in cluster 
galaxies. The SN la rate in spirals drops at z ~ 1.9 because of the low- metallicity 
inhibition of SNe la. We can test the metallicity effect by finding this drop of 
the SN la in spirals, if high-redshift SNe la at z J> 1.5 and their host galaxies are 
observed with the Next Generation Space Telescope. In ellipticals, the chemical 
enrichment takes place so early that the metallicity is large enough to produce SNe 
la at z > 2. The two peaks of SN la rates at z ~ 2.6 and z ~ 1.6 come from 
the MS+WD and the RG+WD systems, respectively. The SN la rate in ellipticals 
decreases at z ~ 2.6, which is determined from the shortest lifetime of SNe la of 
~ 0.5 Gyr. Thus, the total SN la rate decrease at the same redshift as ellipticals, 
i.e., z ~ 2.6. (Note, the decrease of the SN la rate at z ~ 1.6 disappears if we adopt 
Zf ~ 3, because the peak from the MS+WD systems moves to lower redshifts.) 

B In Field 

We also predict the cosmic supernova rates assuming that the formation of el- 
lipticals in field took place for over the wide range of redshifts, which is imprinted 
in the observed spectra of ellipticals in the Hubble Deep Field [64]. The adopted 
SFRs are the same as the case of cluster galaxies, but for the formation epochs Z{ 
of ellipticals distribute as oc exp(— ((z — 2)/2) 2 ) in the range of < z < 5, 

The synthesized cosmic SFR has a broad peak around z ~ 3, which is in good 
agreement with the recent sub-mm observation [65]. Figure 9 shows the cosmic 
supernova rates in field galaxies. As in Figure 8, the SN la rate in spirals drops at 
z ~ 1.9. The averaged SN la rate in ellipticals decreases at z ~ 2.2 as a result of 
~ 0.5 Gyr delay of the decrease in the SFR at z > 3. Then, the total SN la rate 
decreases gradually from z ~ 2 to z ~ 3. 

The rate of SNe II in ellipticals evolves following the SFR without time delay. 
Then, it is possible to observe SNe II in ellipticals around z ~ 1. The difference in 
the SN II and la rates between cluster and field ellipticals reflects the difference in 
the galaxy formation histories in the different environments. 




C Summary 

(1) In the cluster environment, the predicted cosmic supernova rate suggests that 
in ellipticals SNe la can be observed even at high redshifts because the chemical 
enrichment takes place so early that the metallicity is large enough to produce SNe 
la at z > 2.5. In spirals the SN la rate drops at z ~ 2 because of the low-metallicity 
inhibition of SNe la. 

(2) In the field environment, ellipticals are assumed to form at such a wide range 
of redshifts as 1 < z < 4. The SN la rate is expected to be significantly low beyond 
z > 2 because the SN la rate drops at z ~ 2 in spirals and gradually decreases 
from z ~ 2 in ellipticals. 
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